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ABSTRACT 


Radiation  patterns  lnbi  bent  computed  for  an  axial  slot  on 
an  infinite  circular  cylinder  coated  with  a  radially  inhomogeneous 
plasma  sheath.  The  relative  dielectric  constant  was  assumed  to 
have  three  types  of  radial  dependence:  monotonic  with  a  minimum 
at  the  inner  edge  of  the  sheath  but  with  no  portion  of  the  sheath  be¬ 
low  plasma  resonance;  monotonic  with  the  inner  portion  of  the 
sheath  below  plasma  resonance;  parabolic  with  a  minimum  within 
the  sheath  and  the  inner  portion  of  the  sheath  below  resonance. 
The  radial  wave  functions  within  the  sheath  are  obtained  in  the  form 
of  infinite  series.  The  resulting  radiation  patterns  are  found  to  be 
sharply  directed  in  the  forward  direction  relative  to  the  equivalent 
free-space  pattern.  These  sharply  forward-directed  patterns  are 


consistent  with  the  interpretation  of  a>*diffusion-type'®rproces8 


V 


within  layers  of  the  sheath  below  cutoff  and  a  predominantly 
^fcref raction -type  ’^process  within  layers  above  cutoff. 
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I.  INTRODUCTION 


The  recent  literature  contains  numerous  analyses  of  radiation 
from  slot  antennas  on  dielectric -clad  and  plasma-clad  cylinders.  These 
analyses  are  especially  useful  in  predicting  the  requirements  for  com¬ 
municating  with  high-speed,  re-entry  vehicles  of  cylindrical  shape.  With 
few  exceptions  the  sheath  has  been  represented  by  a  homogeneous  coating 
with  an  equivalent  dielectric  constant,  determined  by  the  "average" 
properties  of  the  sheath.  It  is  to  be  expected,  however,  that  important 
effects  caused  by  density  gradients  which  exist  in  re-entry  sheaths  will 
seriously  modify  the  radiation  characteristics  of  the  homogeneous  model. 

Consequently  a  study  of  the  radiation  pattern  of  a  slotted  cylinder 
clad  with  a  radially  inhomogeneous  sheath  has  been  undertaken.  The 
assumed  radial  variation  of  the  dielectric  constant  is  highly  idealized 
(parabolic),  yet  it  corresponds  closely  to  sheath  conditions  of  consider¬ 
able  physical  interest.  It  can  be  shown  that  the  equivalent  dielectric 
constant  of  a  cold  plasma  is  approximately: 

=  1  NeZ 

€0  m  e  Q  u»2 

where  N  is  the  electron  density.  Estimates  of  the  radial  distribution  of 
electron  density  in  the  plasma  sheath  surrounding  a  hypersonic  re-entry 
vehicle  indicate  a  maximum  density  at  or  near  the  surface  of  the  vehicle, 
and  a  monotonically  decreasing  density  with  increasing  distance  from  the 


maximum  (Ref.  1).  Such  a  negative  gradient  of  electron  density  would 
cause  the  dielectric  constant  to  increase  with  increasing  radius  from  a 
value  less  than  unity  to  the  free -space  value.  Furthermore,  it  may  not 
be  unreasonable  to  assume  that  over  a  range  of  radii  the  electron  density 
is  sufficiently  large  that  below  certain  frequencies  the  dielectric  constant 
is  negative  and  that  region  is  opaque  to  propagating  waves.  By  appro¬ 
priate  choice  of  parameters  the  radial  variation  of  dielectric  constant 
assumed  in  the  following  section  can  be  consistent  with  both  relatively 
steep  negative  gradients  of  electron  density  and  portions  of  the  sheath 
below  cutoff. 

The  equivalent  dielectric  constant  used  in  the  following  analyses 
is  derived  with  the  restriction  that  nonlinear  effects,  thermal  effects,  and 
the  effects  of  collisions  may  be  neglected.  It  has  been  pointed  out  by 
Samaddar  in  the  analysis  of  a  similar  problem  (Ref.  2)  and  by  others 
(Refs.  3  and  4)  that  this  approximation  is  no  longer  valid  if  the  dielectric 
constant  becomes  vanishingly  small;  under  such  conditions  it  becomes 
necessary  to  use  the  kinetic  theory  of  plasma  behavior.  Although  several 
"of  the  numerical  examples  considered  allow  the  dielectric  constant  to  pass 
through  zero  at  some  radius  within  the  sheath,  kinetic  theory  has  not  been 
applied.  Consequently  the  interpretation  of  some  of  the  curves  is  limited 
by  the  possibility  that  nonlinear  and  thermal  effects  in  the  vicinity  of 
vanishing  dielectric  constant  may  modify  the  patterns. 
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II.  FORMULATION  OF  THE  PROBLEM 


An  infinite,  perfectly  conducting  circular  cylinder  of  radius  a  is 
coaxial  with  the  z-axis  (Fig.  1).  The  cylinder  is  covered  by  an  inhomo¬ 
geneous,  isotropic  plasma  coating  of  permittivity  «(kp)  =  e0fl  -  — 1 

L  ml0wz  J 

and  permeability  jjl 0  extending  from  a  <  p  <  b.  N(kp)  is  the  electron 
density,  a  function  of  radius.  An  infinite  axial  slot  on  the  surface  of  the 
cylinder  (p  =  a)  extends  from  -  ^  <  $  <  +  Because  of  the  assumed 
nature  of  the  fields  in  the  slot,  the  tangential  electric  fields  at  p  =  a  are: 

Ez(a,0)  =  0  (1) 

'  I0l>^  (2a) 


Ef(a.ff)  = 


(2b) 


Within  the  plasma  coating  from  a  <  p  <  b,  the  permittivity  is 
assumed  to  have  the  following  dependence  on  radius: 

e(p)  =  c0f(kp)  (3) 

where 

f(kp)  =  H(kp  -  kc)  (kp  -  kcA)  =  H(kp  -  kc)^  +  H(kp  -  kc)  (kc  -  kcA)  (4) 
Consequently  f(kp),  the  relative  dielectric  constant,  has  three  types  of 
radial  dependence:  monotonic  with  a  minimum  at  p  =  a  but  with  no  portion 
of  the  sheath  below  plasma  resonance  (Fig.  2a);  monotonic  with  the  inner 
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portion  of  the  sheath  below  plasma  resonance  (Fig.  2b);  parabolic  with 
a  minimum  at  pQ  >  a  and  the  inner  portion  of  the  sheath  below  plasma 
resonance  (Fig.  2c).  In  each  case  the  outer  edge  of  the  sheath  is  defined 
by  € (b)  =  C0. 

The  corresponding  radial  variations  of  electron  density  are  also 
indicated  in  Figs.  2a,  2b  and  2c. 


3 


s  Fig.  1  Cross  Section  Of  Infinite  Slotted 

Cylinder  With  Plasma  Sheath. 

I 
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III.  SOLUTION  OF  MAXWELL'S  EQUATIONS 


J  IN  THE  INHOMOGENEOUS  MEDIUM 


i 


> 


-iwt 


For  an  assumed  harmonic  time  dependence  e  (suppressed 
throughout),  Maxwell's  equations  become: 


V  x  H  =  =  -iuf(p)  E 

0  t 


(5a) 


Vx  E 


a_B 

dt 


+  iw 


^o 


H 


(5b) 


The  electric  field  may  be  eliminated,  yielding 

yx  y  x  H  =  x  (Vx  H)  +  u2pi0  «Qf  (kp)H  (6) 

Since  no  radial  component  of  the  magnetic  intensity  will  be  excited,  it  is 
expected  that  H  can  be  expressed  as  the  curl  of  a  radially  directed 
vector  function  of  p  and  ft  (after  the  vector  wave -function  method  of 
Stratton  and  Hansen  [Refs.  5  and  6]): 

H  =  Vx  {(kp)  R(kp)  H(0)  Ip}  (7) 

Substitution  of  this  quantity  into  Eq.  6  yields 

(8) 

Equation  8  is  satisfied  by  requiring  the  bracketed  quantity  itself  to  be 


=0 
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zero.  This  condition  then  yields  the  two  equations  defining  R(kp)  and  H(0): 

H"  (0)  +  m2  H(0)  =  0  (9) 

{(kp)2}R‘'(k/9)+  } R‘ (kpH- {(kpfflkp)- m*}  R  ( k^> )  =0  (10) 

The  solutions  of  Eq.  9  are  clearly 

H(jf)  =  Cm  cosm^  +  Dj^sinm^  (11) 

The  solutions  of  Eq.  10  are  two  independent  radial  functions  which  will  be 
determined  in  a  later  section: 

R(kp)  =  Am  R^kp)  +  Bm  R^(kp>  (12) 

In  free  space  f(kp)  =  1  and  Eq.  10  becomes  Bessel's  equation 

(kp)2  R"  (kp)  +  (kp)  R'  (kp)  +  f  (kp)2  -  m2J  R(kp)  =  0  (10a) 

In  this  special  case  RmW)  and  R^(kp)  are  the  standard  Bessel  or 
Hankel  functions  of  integral  order. 


IV.  MATCHING  THE  BOUNDARY  CONDITIONS 


For  the  geometry  of  the  plasma-clad  axially  slotted  cylinder 
*  described  previously,  the  fields  may  be  expanded  in  the  following 

\,  manner: 
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Et(p’f)=zr! h  m^0  ^[AmR™)(kp)+ BmR“)(kp)]C0,m^  (13a) 

*<p£  b 

Hz(p,*)  =  -k  X  [An.Rii^kp)  +  BrnR^)(kp)]cosm^  (13b) 

m=o  L  J  ' 

00  1 
Eflp>  A  =  377  m?0  ^  [Cm  Hm)(kf)] co8m* 


CD 


HZ(P,0)  =  -k  £  CmH^(kp)  cosm^ 
m  =0 


p  >b 


(14a) 


(14b) 


Matching  the  boundary  conditions  on  the  tangential  field  at  p  =  a  and 
p  =  b  and  letting  A  — »0,  it  can  be  shown  that 


'm 


"i  '  ^m3] 

(*+8om)  (^m2)  +  j 


21  where  Gq  = 


w  c(a)  Vo 
it  ak^ 


(15) 


and 


*ml  =  Rm'  (kb)  R-m>(kb)  -  R^)’  (kb)  R^  (kb)  (16) 

‘ml  -  xi/  (ka|{jm(kl.)[^  R^)(kb)-R^)'  (kb)]  -  Jm+1(kb)R^’(kb|} 

-  B^1’  (ka){jm(kb>[^  <>(kb)  -  R^1'  (kb)]  -  Jmtl(kb)  R^'(kb)]  (17a) 
'm3  *  *U?  (k»{Nm(kb)[gR^)(kb)  -  R^’  (kb)]  -  Nmtl(kb)  R^’(kb)} 

-Bm’1  (kal{Nm(kb)|]^  R^’tkb)-  R^1  (kb)]  -  Nm+1(kb)  R^’lkb)}  (17b) 
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4 


Upon  substitution  of  Eq.  15  and  the  asymptotic  value  of  the  Hankel 


i 

/ 


function  into  Eq.  14b,  the  far-zone  field  becomes: 


Hz(kp) 


3ir  _ 
Mkp  -  T)  ® 

e  L 

m=o 


mir 

e  fml(fm2-ifm3>  cogmjf 


(18) 


V.  SOLUTION  FOR  THE  RADIAL  FUNCTIONS 


i 

i 


\ 


Inspection  of  Eq.  10  reveals  that  p  =  c  is  a  regular  point  of  the 
differential  equation.  It  is  then  to  be  expected  that  R^(kp)  and  R^(kp) 
can  be  expanded  as  an  infinite  power  series  about  p  =  c  (Ref.  7). 


Equation  10  can  then  be  rewritten: 

A  (kp)  R"  (kp)  +  B (kp)  R'  (kp)  +  C  (kp)  R(kp)  =  0  (19) 

where 

* 

A  (kp)  =  Y  an  (kp  -  kc)n  (20) 

n  =  l 

aj  =  H(kc)2  (kc  -  kcA>  (20a) 

a^  =  H(kc)  (3kc  -  2kc^)  (20b) 

a3  =  H(3kc  -  kc^)  (20c) 

a4  =  H  (20d) 
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and 


s 

=  £  b^kp  - kc)n  (21) 

n  =  o 

bo  =  -H(kc)2(kc  -  kcA)  =  -a!  (21a) 

bi  =  -H(kc)  (3kc  -kcA)  (21b) 

b2  =  -3H  (kc)  (21c) 

b3  =  -H  =  -a4  (21d) 

and 

« 

C(kp)  =  Y  cn  (kp  *  kc  )n  (22) 

ns  i 

ci  =  -m2H(kc-kcA)  (22a) 

c2  =  H^H(kc)2(kc -kcA)2  -  m2J  (22b) 

c3  =  2H2  (kc)  (2kc  -  kcA)  (kc  -  kcA)  (22c) 

c4  =  H2^6(kc)2  -  6(kc)(kcA)+  (kcA)2J  (22d) 

c5  =  2H2(2kc-kcA)  (22e) 

c6  =  H2  (22f) 

If  Rm(kp)  is  assumed  to  be  of  the  form 
<9 

Rm(kp)  =  £  dn(kp-kc)n+a  (23) 

n=o 

then  substitution  into  Eq.  19  yields  the  following  indicial  equation: 

a[ax  (a-1)  +  bQJ  =0  (24) 
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The  solutions  of  Eq.  24  are 


> 


t 

( 


\ 


a  -  0 


(24a) 


(24b) 


Following  the  procedure  outlined  in  Ref.  7  if  a  solution  is  assumed  using 
the  algebraic  smaller  of  the  two  indices,  two  independent  infinite  series 
are  generated  if  m  =  0.  It  can  be  shown  that  if  m  =  0: 

R^Nkp)  =  X  -kc)n+2  (25) 

n  =  o 


where 

0d01)sl  (25a) 

0  dn  '  =  -{(n+l)(na2+b1)()  d^i  +  n  [(n-1)  a3  +  b2]  Q  d^.2 +  [(n-l)(n-3)  a4+  c2J  Q  dj^ 


+  c  - 


0dn 


(1)  Jl) 

4+  c4  0  dn-5  + 


c5 


rf(1) 
0  dn-6 


+  c6  0  dn -?}/(*»+ 2)n  a  j,  n  >  1  (25b) 


and 


RqVp)  = 


X  0dn2^(kp-kc)n 
n=o 


where 


1 


0 


0 


(26) 


(26a) 

(26b) 

(26c) 
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(26d) 


t 


0dn2)  =  -{(n-l)[(9-2)a2+b1J0di2.)i  +(n+2)[tn-3)a3  +  b2J0dJl2.)2 

r  i  (2)  it)  J2) 

+  |(n-3)  (n-5)a4  +  c2JQdn_3  +  c3  0dn_4  +  c4  0dn.5 

+  c5  0di2-6  +  c6  0d  k2-7}/n<n-2>al  *  n  >  3 

In  the  case  m  ^  0,  both  roots  of  the  indie ial  equation  yield  the 


same  series: 


Rm  (kp)  =  £  mdn  '(kp-kc) 

n=  o 


CD 


.(1). 


where 


.(1) 


n+2 


(27) 


m  a0  =  1 


(27a) 


m ' 


H 

^(n+l)(na2+b1)  m 

4"  ♦[> 

i[jn- l)a3  +  bjf  c ; 

1  H(1) 

Ljmdn-2 

(n-l)(n-3)a4+  c2 

,  id 

(1) 

(1) 

+ 

jrndn-3+c3rridn-4+c4m 

dn-5+  c5 

U) 


+  C6md^-7}/n(n+<i)al  ’  n-1 


(1) 


(27b) 


Rq  ^  (kp)  is  the  limiting  case  of  R^(kp)  with  m  =  0.  A  second  inde¬ 
pendent  solution  can  be  generated  of  the  form  (Ref.  7) 

00 


R^(kp)  =  R^(kp)  -Mkp -kc)Z  +  £  mfn(kp-kc)n 


(28) 


n  =  o 


Substitution  of  Eq.  28  into  Eq.  19  then  yields  the  recursion  relations  for 


m  *n: 


m 


f0  *  *{4al/cl} 


(28a) 
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0 


m*  i  = 


(28b) 


mL2 


f,  =  0 


(28c) 

(U 


m 


fn  =  - {4(n-  l)ai  m  £}2  +  2[(2n-3)a2  +  bxj  m  +  2  [(2n-5)a3  +  b2]  m d'nV4 

+  4(n-4)a4md|11.5+  (n- l)£(n-2)a2  +  bxj  mfn_  x  +  |jn-2)[(n- 3)a3  +  b2]+ c  xJm  fn_2 
+  [(n-3)(n-5)a4+  c2]m*n-3+c3  m^n-4  +  c4m*n-5  + c5m^n-6 


+  c6mfn-7}/n<n-2>al* 


(28d) 


Equations  25,  26,  27,  and  28  can  then  be  substituted  into  Eq.  18 
to  yield  the  far-field  pattern.  Considerable  simplification  can  be  accom¬ 
plished  by  substitution  of  the  Wronskian  relations,  developed  in  Appendix  A, 
into  the  formulas.  The  far-field  pattern  then  becomes 


P(0)  %/[R(*>]2  +[l  <^]2 


(29) 


where 


nrt IMM-MWt-Mtlf  Jci. _ +e(kc)t?'  g-S^emgr'mJ  coamtf> 

(kblte-kc.) 


(30) 


1«*»  = 


(kc)(kb-kcXkb 

(kbXkc-kcA) 


-kc4)f  f os  nf.  tmg+coi^fm,] 

— T7 — .8.  8(kc)  >  k - - - = - J—  cotmcp 

I  l  <WW  9 


(3D 
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In  the  case  of  a  vanishing  sheath,  i.  e. ,  a  =  b,  it  can  easily  be 
shown  that  Eq.  29  reduces  to 

.  .  oo  -im  ~  .  , 

P(0)  =  £  e  1  cQsm^  (32) 

m'°  d+W^'tkb) 

which  is  the  well-known  result  for  a  slotted  cylinder  in  free  space. 


VI.  NUMERICAL  EVALUATION  OF  THE  FIELDS 


i 


\ 


The  magnitude  of  the  far-zone  field  pattern  (cf.  Eq.  18)  has  been 
evaluated  on  a  large  digital  computer  for  various  configurations  of 
dielectric  constant.  These  patterns  are  plotted  in  the  following  figures. 

An  inner  sheath  radius  of  ka  =  5  and  an  outer  radius  of  kb  =  6  have  been 
.maintained  throughout.  Only  half  of  each  pattern  has  been  plotted  since 
symmetry  exists  about  0  =  0.  Each  pattern  has  been  normalized  to  its 
maximum  value  (consistently  in  the  forward  direction).  Consequently, 
direct  amplitude  comparisons  may  not  be  made  between  individual  patterns. 
However,  comparison  of  relative  pattern  shapes  has  been  facilitated  by 
this  normalization. 

The  field  patterns  plotted  in  Fig.  3  resulted  from  a  dielectric 
constant  variation  of  the  type  indicated  in  Fig.  2a:  a  monotonic  variation 
with  a  minimum  at  the  inner  edge  of  the  sheath  but  with  no  portion  of  the 
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Fig.  3  Field  Patterns  For  Slotted  Cylinder  With  Inhomogeneous  Sheath  Above  Resonance. 

a.  ka  =  5.0.  kb  =  6.0,  kc  =  4.0,  kcA  =  1.5,  H=  0.1111 

b.  ka  =  5.0,  kb  =  6.0,  kc  =  4.5,  kcA  =  2.5,  H  =  0.1905 

c.  ka  =  5.0,  kb  =  6.0,  kc  =  4.8,  kcA  =  3.0,  H  =  0.2778 


sheath  below  plasma  resonance.  The  corresponding  radial  dependence 
of  dielectric  constant  and  associated  electron  density  are  plotted  in  Fig.  4. 
The  free-space  pattern  of  a  slotted  cylinder  without  any  sheath  has  been 
included  for  reference.  This  pattern  exhibits  ripples  at  large  angles 
which  may  be  interpreted  as  interference  among  circumferentially 
directed  waves.  When  an  inhomogeneous  sheath  is  present  these  ripples 
no  longer  exist,  with  the  exception  of  a  small  local  maximum  at  180  deg. , 
a  characteristic  feature  of  many  types  of  diffraction  patterns.  In  addition, 
the  sheath  patterns  indicate  increased  power  radiated  in  forward  directions, 
relative  to  the  free-space  pattern,  with  a  correspondingly  reduced  power 
radiated  toward  the  rear.  It  is  evident  from  a  comparison  of  the  three 
patterns  that  greater  forward  enhancement  results  from  steep  sheath 
gradients. 

Although  within  the  sheath  the  condition 


— —  ilM  «  i 

k€(p)  dp 


(33) 


is  not  satisfied,  the  resulting  field  patterns  in  Fig.  3  appear  to  be 
consistent  with  a  quasi-geometrical  ray- optics  interpretation.  Inasmuch 
as  the  dielectric  constant  increases  with  increasing  radius,  each  "ray" 
emerging  from  the  slot  is  refracted  continuously  in  the  forward  direction, 
thus  reducing  or  eliminating  the  effects  of  circumferential  waves  and 
producing  field  patterns  enhanced  in  the  forward  directions.  The  steeper 


sheath  gradients  produce  greater  bending  of  the  "rays"  in  the  forward 
direction  and,  consequently,  greater  forward  enhancement.  Although 
this  interpretation  appears  valid,  the  steep  gradients  within  the  sheaths 
undoubtedly  produce  other  effects  (e.g.  scattering)  in  addition  to 
refraction,  all  of  which  are  included  in  Eq.  18,  the  complete  solution  of 
Maxwell's  equations. 

The  field  patterns  plotted  in  Figs.  5a,  6a,  and  7a  have  been 
derived  for  inhomogeneous  sheaths  that  are  partially  below  plasma 
resonance.  The  corresponding  dielectric  constant  variations  are  plotted 
in  Figs.  5b,  6b,  and  7b,  respectively.  In  each  case  the  minimum  value  of 
dielectric  constant,  with  the  associated  maximum  electron  density,  occurs 
at  the  inner  edge  of  the  sheath,  p  =  a.  It  is  evident  that  (1)  wide-angle 
ripples  no  longer  exist,  and  (2)  these  patterns  are  more  sharply  forward- 
directed  than  either  the  free-space  pattern  or  the  patterns  considered 
previously  for  inhomogeneous  sheaths  above  plasma  resonance  (cf.  Fig.  3). 

Propagating  waves  cannot  exist  in  the  region  below  plasma 
resonance.  Instead  the  fields  are  evanescent  in  nature,  decaying  strongly 
along  each  "ray  path".  Consequently  the  strongest  "rays"  emerging  from 
the  cutoff  region  will  have  undergone  the  shortest  "ray  paths",  resulting  in 
a  sharply  forward-directed  pattern.  This  effect  has  been  previously  calcu¬ 
lated  for  homogeneous  sheaths  below  plasma  resonance  (Ref.  8).  The  "rays" 
will  undergo  additional  forward  bending  upon  emerging  into  the  outer  region 
that,  although  above  cutoff,  is  radially  inhomogeneous.  Consequently,  it  is 
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Fig.  5b  Radial  Dependence  Of  Dielectric  Constant  And  Electron 
Density  For  Inhomogeneous  Sheath  Partially  Below 

Resonance. 


ka  =  5.0,  kb  =  6.0,  kc  =  5.3,  kc^  =  4.3,  H  =  0.8403 
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Fig.  6a  Field  Pattern  For  Slotted  Cylinder  With  Inhomogeneous  Sheath  Partially  Below  Resonance, 
ka  =  5.0,  kb  =  6.0,  kc  -  5.5,  kcA  =  4.5,  H  =  1.3333 
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Fig.  6b  Radial  Dependence  Of  Dielectric  Constant  And  Electron 
Density  For  Inhomogeneous  Sheath  Partially  Below 

Resonance. 


ka  =  5.0,  kb  =  6.0,  kc  =  5.5,  kc^  =  4. 5,  H  =  1.3333 
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Fig.  7a  Field  Pattern  For  Slotted  Cylinder  With  Inhomogeneous  Sheath  Partially  Below  Resonance, 
ka  =  5.0,  kb  =  6.0,  kc  =  5.8,  kc^  =  4.2,  H  =  2.7778 
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Fig.  7b  Radial  Dependence  Of  Dielectric  Constant  And  Electron 
Density  For  Inhomogeneous  Sheath  Partially  Below 

Resonance. 

ka  =  5.0,  kb  =  6.0,  kc  =  5.8,  kcA  *  4.2,  H  *  2.7778 
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not  unreasonable  to  assume  that  both  the  "diffusion- type 11  process  within 
the  inner  region  and  the  " refraction- type  "  process  within  the  outer  will 
contribute  to  a  relatively  sharp  maximum  in  the  forward  direction. 
Comparison  of  Figs.  5-7  with  Fig.  3  indicate  that  processes  within  the 
cutoff  region  may  contribute  more  to  the  high  degree  of  forward  direc¬ 
tivity  than  processes  within  the  outer. 

However,  the  electromagnetic  effects  taking  place  within  inhomo¬ 
geneous  sheaths  of  this  type  are  extremely  complex.  No  apparent 
correlation  is  evident  between  pattern  features  and  sheath  distributions 
for  the  three  examples  considered  in  Figs.  5-7.  Of  the  three,  Fig.  6a  is 
the  most  directive,  although  the  corresponding  sheath  for  this  pattern  is 
intermediate  in  both  cutoff  layer  thickness  and  sheath  gradient.  Since  the 
dielectric  constant  becomes  zero  within  the  sheaths  considered  in  Figs.  5-7, 
it  is  evident  from  Eq.  33  that  the  approximations  of  geometrical  ray  optics 
are  even  less  valid  than  previously.  Furthermore,  increased  gradients 
may  result  in  more  significant  scattering  effects.  A  simple  interpretation 
is  not  readily  apparent  for  sheaths  of  this  type. 

The  field  pattern  in  Fig.  8a  was  also  derived  for  an  inhomogeneous 
sheath  partially  below  plasma  resonance.  However,  the  minimum  value 
of  dielectric  constant  occurs  within  the  cutoff  region  at  pQ  where  pQ  >  a, 
as  indicated  in  Fig.  8b.  Again,  there  is  a  high  degree  of  directivity  in  the 
forward  direction.  However,  there  is  no  significant  difference  between 
this  pattern  and  those  considered  in  Figs.  5-7. 
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Fig.  6a  Field  Pattern  For  Slotted  Cylinder  With  Inhomogeneous  Sheath  Partially  Below  Resonance 
ka  =  5.0,  kb  =  6.0,  kc  =  5.5,  kc^  =  4.8,  H  =  1.6667 
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Fig.  8b  Radial  Dependence  Of  Dielectric  Constant  And  Electron 
Density  For  Inhomogeneous  Sheath  Partially  Below 
Resonance. 


ka  =  5.0,  kb  =  6.0,  kc  =  5.5,  kc^  -  4.8,  H  =  1.6667 
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VII.  SERIES  CONVERGENCE 


The  circle  of  convergence  for  each  of  the  radial  series  (Eqs.25, 

26,  27,  28)  is  defined  by  |kp  -  kc  |<(kc  -  kc^)  (Ref.  7).  These  series 
converge  extremely  rapidly  if  j  kp  -  kc  |<  I.  This  condition  was  satisfied 
in  the  computation  of  Figs.  5-8,  i.  e. ,  when  a  portion  of  the  sheath  was 
below  plasma  resonance.  It  was  found  that  for  m  ^  34  (the  greatest 
value  of  m  considered)  the  35-term  and  99-term  series  were  indistin¬ 
guishable  (to  the  four  significant  figures  of  the  computer  printout  roundoff). 
As  an  additional  check,  the  Wronskian  relation  agreed  to  within  four 
significant  figures.  In  each  case  the  field  patterns  were  identical  for  the 
35-term  and  99-term  series. 

In  cases  where  |  kp  -  kc |  ^  1,  a  discrepancy  in  the  fourth  significant 
figure  was  noticed  between  the  35-term  and  99-term  series  for  m  >  12. 
This  resulted  in  a  slight  variation  (in  the  fourth  significant  figure)  in  the 
small  values  of  the  field  patterns  between  the  35-term  and  99-term  results. 
However,  the  field  patterns  were  identical  for  the  60-term  and  99-term 
series,  and  the  Wronskian  relations  again  agreed  within  the  printout 
roundoff  for  the  35-term  and  99-term  series.  It  was  evident  that  by  taking 
60  terms  or  more  all  series  converged  to  four  significant  figures  for  the 
range  of  parameters  considered.  In  cases  where  thicker  sheaths  or  sheaths 
of  larger  radius  are  considered,  convergence  may  be  a  matter  of  increased 
significance. 
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APPENDIX  A 


DERIVATION  OF  THE  WRONSKIAN  RELATIONS 


Since  both  R^(kp)  and  R^(kp)  satisfy  Eq.  10,  the  following  pair 


of  equations  can  be  formed: 

{(|V>>2}Rfn(lS0>+-  {f(tyi) -(kyo)  }Rm(iy>)  +  j(k/>f  %>) - rr? }  0y>)  =0 

(Al) 

{(K/of}Rm<VJ"l_  1T^5  }»WV),+{(,y>)  Rm  « kP>  =0 

(A2) 

Multiplying  (Al)  by  R^(kp)  and  (A2)  by  R^(kp)  and  subtracting 

yields: 

frpf  [r^ (*Vo)R'Si (k/>)-  Oyo)  Rm1*4  (kyo)]+  [ftlpHkjo)^-^] 

[  Rm  (k p)  R?f  (k/> )  -  R„ 1  (l ip)  Rm<2)  ( k/> )  j  »  0 

(A3) 

This  can  be  rewritten  as 


d 

d(kp) 


in  {r^1  (kp)  R^'lkp)  -  (kp)  r|„2)  '(kp)}  =  ^  tn{®*]}  (A4) 


(kp) 


Integrating: 


Af(kp) 

(kp) 


(A5) 


The  constant  of  integration,  A,  can  then  be  evaluated  by  taking  the  leading 
terms  of  the  left-hand  side  of  Eq,  A5  as  p  -»  c.  Hence 


m  =  0: 


m  £  1: 


A  (2kc) 

A  '  H(kc-kcA) 

A  -  8(kc)3 

m^  H(kc  -  kc^) 


(A6) 

(A7) 
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